Dynamical control of matter-wave tunneling in periodic potentials 
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We report on measurements of dynamical suppression of inter-well tunneling of a Bose-Einstein 
condensate (BEC) in a strongly driven optical lattice. The strong driving is a sinusoidal shaking 
of the lattice corresponding to a time-varying linear potential, and the tunneling is measured by 
letting the BEC freely expand in the lattice. The measured tunneling rate is reduced and, for certain 
values of the shaking parameter, completely suppressed. Our results are in excellent agreement with 
theoretical predictions. Furthermore, we have verified that in general the strong shaking does not 
destroy the phase coherence of the BEC, opening up the possibility of realizing quantum phase 
transitions by using the shaking strength as the control parameter. 

PACS numbers: 03.65.Xp, 03.75.Lm 
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Quantum tunneling of particles between potential wells 
connected by a barrier is a fundamental physical effect. 
While typically quantum systems decay faster when they 
are perturbed, if the wells are periodically shaken back 
and forth (or a time- varying potential is applied in a dif- 
ferent way), the tunneling rate can actually be reduced 
and, for certain shaking strengths, even completely sup- 
pressed [HQ- 

Modifications of the dynamics of quantum systems by 
applying periodic potentials have been investigated in 
a number of contexts including the renormalization of 
Lande g-factors in atoms Q , the micromotion of a single 
trapped ion Q and the motion of electrons in semicon- 
ductor superlattices [1| . In particular, theoretical studies 
of double- well systems and of periodic potentials have led 
to the closely related concepts of coherent destruction of 
tunneling and dynamical localization [JQ. In the latter, 
tunneling between the sites of a periodic array is inhib- 
ited by applying a periodically varying potential, e.g. by 
shaking the array back and forth (see Fig. 1), and as a 
consequence the tunneling parameter J representing the 
gain in kinetic energy in a tunneling event is replaced by 
\Jcs\ < l"^!- In a number of experiments signatures of 
this tunneling suppression have been observed d, 0, [1] , 
and recently dynamical localization and coherent sup- 
pression of tunneling have been demonstrated using light 
propagating in coupled waveguide arrays Also, 
the predictions of the Bose-Hubbard model in a moving 
frame were recently tested So far, however, an ex- 



act experimental realization of the intrinsically nonlinear 
Bose-Hubbard model Q driven by a time-periodic poten- 
tial has not been reported. 

In this Letter, we report on the observation of the 
dynamical tunneling suppression predicted in refs. [2, 
Il2i] using Bose-Einstein condensates (BECs) in strongly 
driven periodic optical potentials In contrast to 

other systems, the characteristics of such optical lattices 
- potential depth, lattice spacing, driving strength and 
frequency - can be freely chosen and allow us to control 
the tunneling over a wide range of parameters. In this 
way we were able to experimentally confirm theoretical 
predictions with great accuracy. Also, our system allows 
us to observe the effects of the shaking both by monitor- 





FIG. 1: Suppression of tunneling by strong driving. The 
dynamics of a Bose-Einstein condensate in a periodic po- 
tential is governed by the tunneling matrix element J and 
the on-site interaction energy U (above). If the potential is 
strongly shaken, tunneling between the wells is dynamically 
suppressed, leading to a renormalized tunneling matrix ele- 
ment Jeff (below) but leaving the interaction energy U unaf- 
fected. 



ing the real-space expansion of the BEC in the optical 
lattice and by performing time-of-flight experiments in 
which the phase coherence of the BEC can be measured. 
The latter experiments allow us to verify that the tunnel- 
ing suppression occurs in a phase-coherent way in spite 
of the strong shaking. 

Furthermore, BECs have an intrinsic nonlinear on-site 
interaction energy (represented by U in Fig. 1), the in- 
terplay of which with the tunneling parameter J has been 
shown to lead to the Mott-insulator quantum pha se tran- 
sition for a critical value of the ratio U/J [13 , hj. It has 
been theoretically predicted that for a BEC in a shaken 
optical lattice, this ratio can be replaced by U / JcS and 
hence that it should be possible to drive the system across 
the quantuin phase transition by varying the shaking pa- 
rameter [2, . In this work, we demonstrate the feasi- 
bility of the key ingredients of this scheme. In particu- 
lar, we show that when tunneling in the shaken lattice is 
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completely suppressed, the phase coherence of the BEC 
is lost in agreement with the physical picture of a sudden 
'switch-off' of the inter-well coupling and a subsequent 
independent evolution of the local phases due to colli- 
sions between the atoms [H, |T3] ■ 

Our system consisting of a Bose-Einstcin condensate 
inside a (sinusoidally) shaken one-dimensional optical lat- 
tice is approximately described by the Hamiltonian 

ffo = -J^(c|c,+4cO+|l]"j(^ij-l)+^cos(a;i)^ 

(ij) j j 

(1) 

where c'p are the boson creation and annihilation opera- 
tors on site fii = c\ci are the number operators, and K 
and Lo are the strength and angular frequency of the shak- 
ing, respectively. The first two terms in the Hamiltonian 
describe the Bose-Hubbard model fli'l with the tunnel- 
ing matrix element J and the on-site interaction term 
U . The shaking of the lattice is expected to renormalize 
the tunneling matrix element J, leading to an effective 
tunneling parameter 



time- varying force 

F{t) = mwdiAz^niax cos(cji) = F^ax cos{ujt). (3) 

The peak shaking force F,„ax is related to the shaking 
strength K appearing in Eq. (1) by 



and hence the dimensionless shaking parameter 



(2) 



where is the zeroth-order ordinary Bessel function and 
we have introduced the dimensionless parameter i^o = 
K/huj. 

In our experiment we created BECs of about 5 x 10'' 
87-rubidium atoms using a hybrid approach in which 
evaporative cooling was initially effected in a magnetic 
time-orbiting potential (TOP) trap and subsequently in 
a crossed dipole trap. The dipole trap was realized using 
two intersecting gaussian laser beams at 1030 nm wave- 
length and a power of around 1 W per beam focused 
to waists of 50 /xm. After obtaining pure condensates 
of around 5 x 10"* atoms the powers of the trap beams 
were adjusted in order to obtain elongated condensates 
with the desired trap frequencies 20 Hz in the lon- 
gitudinal direction and 80 Hz radially). Along the axis 
of one of the dipole trap beams a one-dimensional opti- 
cal lattice potential was then added by ramping up the 
power of the lattice beams in 50 ms (the ramping time 
being chosen such as to avoid excitations of the BEC). 
The optical lattices used in our experiments were cre- 
ated using two counter-propagating gaussian laser beams 
(A = 852 nm) with 120 /im waist and a resulting optical 
lattice spacing = A/2 = 0.426 /xm. The depth Vb of 
the resulting periodic potential is measured in units of 
E'rec = tt'^ / {2md'j^) , where m is the mass of the Rb 
atoms. By introducing a frequency difference Av be- 
tween the two lattice beams (using acousto-optic mod- 
ulators which also control the power of the beams), the 
optical lattice could be moved at a velocity v — dLAv 
or accelerated with an acceleration a — rfi^^. In or- 
der to periodically shake the lattice, Av was sinusoidally 
modulated with angular frequency w, leading to a time- 
varying velocity v{t) = diAz^max sin(u;t) and hence to a 
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The spatial shaking amplitude Axmax can then be writ- 
ten as 



AXr, 



2 LOr, 



-Kodi 



(6) 



so for a typical shaking frequency ujl-K — 3 kHz we have 
Asmax « 0.5dL at ifo = 2.4. 

After loading the BECs into the optical lattice, the 
frequency modulation of one of the lattice beams cre- 
ating the shaking was switched on either suddenly or 
using a linear ramp with a timescale of a few millisec- 
onds. Finally, in order to measure the effective tunneling 
rate | Jcff| between the lattice wells (where the modulus 
indicates that we are not sensitive to the sign of J, in 
contrast to the time-of-flight experiments described be- 
low), we then switched off the dipole trap beam that 
confined the BEC along the direction of the optical lat- 
tice, leaving only the radially confining beam switched 
on (the trap frequency of that beam along the lattice di- 
rection was on the order of a few Hz and hence negligible 
on the timescales of our expansion experiments, which 
were typically less than 200 ms) . The BEC was now free 
to expand along the lattice direction through inter-well 
tunneling and its in-situ width was measured using a res- 
onant flash, the shadow cast by which was imaged onto 
a CCD chip. The observed density distribution was then 
fitted with one or two gaussians. 

In a preliminary experiment without shaking (A'o = 0), 
we verified that for our expansion times the growth in 
the condensate width along the lattice direction was 
to a good approximation linear and that the dependence 
of d(Jx/dt on the lattice depth (up to Vo/-Brec = 9) fol- 
lowed the expression for J(Vq/£'i.oc) in the lowest energy 
band [H 
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This enabled us to confirm that da^/dt measured at a 
fixed time was directly related to J and, in a shaken 
lattice, to \Jcff{Ko)\. The results of our measurements of 
\Jcs{Ko)/J\, for various lattice depths Vq and driving fre- 
quencies Lo are summarized in Fig. 2. We found a univer- 
sal behaviour of \Jcff/J\ that is in very good agreement 
with the Bessel- function re-scaling of Eq. ^ . We were 
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able to measure \Jcs/J\ for Kq up to 12, albeit agree- 
ment with theory beyond Kq « 6 was less good, with the 
experimental values lying consistently below the theoret- 
ical curve. For the zeroes of the j7o Bessel function at 
« 2.4 and 5.4, complete suppression of tunneling was 
observed (within our experimental resolution, we could 
measure a suppression by at least a factor of 25). 
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FIG. 2: Dynamical suppression of tunneling in an optical 
lattice. Shown here are the values of |Jcff/J| as calculated 
from the expansion velocities as a function of the shaking 
parameter Kq. The lattice depths and shaking frequencies 
were: Vb/Sroc = 6, aj/27r = IkHz (squares), Vo/Eicc ~ 6, 
Lu/2-K = 0.5 kHz (circles), and Vo/Erec = 4, uj/2tt = IkHz 
(triangles). The dashed line is the theoretical prediction. In- 
sert: Dependence of the tunneling suppression | Jeff/J| on the 
shaking frequency ui for Ko ~ 2.0 and VoZ-Brcc = 9 corre- 
sponding to J/h — 90 Hz. 

We also checked the behaviour of \Jcs/J\ as a func- 
tion of u for a fixed value of Kq = 2 (see insert in 
Fig. 2) and found that over a wide range of frequen- 
cies between Hlo/ J « 0.3 and Tilo/ J 30 the tunneling 
suppression due to the shaking of the lattice works, al- 
though for fiu)/J < 1 we found that \Jcf{{Ko)/J\ as a 
function of Kq deviated from the Bessel function near 
the zero points, where the suppression was less efficient 
than expected. In the limit of large shaking frequencies 
{lj/2h ^ 3 kHz, to be compared with the typical mean 
separation of ~ 15 kHz between the two lowest two en- 
ergy bands at Vo/i?rcc — 9), we observed excitations of 
the condensate to the first excited band of the lattice. 
In our in-situ expansion measurements, these band exci- 
tations were visible in the condensate profile as a broad 
gaussian pedestal below the near-gaussian profile of the 
ground-state condensate atoms. From the widths of those 
pedestals we inferred that \Jcs/ J\ of the atoms in the ex- 



cited band also followed the Bessel-function rescaling of 
Eq. ([2]), and that the ratios of the tunneling rates in the 
two bands agreed with theoretical models. 
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FIG. 3; Phase coherence in a shaken lattice, (a) Dephas- 
ing time Tdcph (decay time of the visibility) of the conden- 
sate as a function of Ko for a lattice with Vo/Etbc = 9 and 
uj/2-K — 3 kHz. The vertical dashed line marks the position of 
Ko — 2.4 dividing the regions with Jeff > (left) and Jeff < 
(right). In both regions, a typical (vertically integrated) in- 
terference pattern of a time-of-flight experiment without final 
acceleration to the zone edge is shown (on the a;-axis, the 
spatial position has been converted into the corresponding 
momentum in units of the recoil momentum prec ~ h/dL-) 
Insert: Rephasing time after dephasing at Ko = 2.4 and sub- 
sequent reduction of Ko- (b) Dephasing time as a function of 
the normalized driving frequency fiuj/ J for Ko = 2.2. 

We now turn to the phase coherence of the BEC in the 
shaken lattice. In order to quantify the degree of phase 
coherence, after shaking the condensate in the lattice for 
a fixed time between 1 and ~ 200 ms we accelerated the 
lattice for « 1 ms so that at the end of the accelera- 
tion the BEC was in a staggered state at the edge of the 
Brillouin zone. After switching off the dipole trap and 
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lattice beams and letting the BEC fall under gravity for 
20 ms, this resulted in an interference pattern featuring 
two peaks of roughly equal height [l^ . In the region be- 
tween the first two zeroes of the Bessel function, where 
Jq < 0, we found an interference pattern (see Fig. 3 (a)) 
that was shifted by half a Brillouin zone, in agreement 
with theoretical predictions. We then measured the visi- 
bility V = (/lmax-/lmin)/(/j max "I" ^min) of the interference 
pattern as a function of the time the condensate spent 
inside the shaken lattice, where /imax is the mean value 
of the condensate density at the position of the two in- 
terference peaks and /imin is the condensate density in a 
region of width equal to about 1/4 of the peak separa- 
tion centered about the halfway point between the two 
peaks. For a perfectly phase-coherent condensate V ~ 1, 
whereas for a strongly dephased condensate V ~ 0. For 
Kq ^2.2, the BEC phase coherence was maintained for 
several tens of milliseconds, demonstrating that the tun- 
neling could be suppressed by a factor of up to 10 over 
hundreds of shaking cycles without significantly disturb- 
ing the BEC. This result is expressed more quantitatively 
in Fig. 3 (a). Here, the condensate was held in the 
lattice {Vq/Ej-cc = 9), and the shaking was switched on 
suddenly at t = (we found no significantly different 
behaviour when Kq was linearly ramped in a few mil- 
liseconds). Thereafter, the visibility was measured as a 
function of time and the decay time constant Tdcph of the 
resulting near-exponential function was extracted. Apart 
from a slow overall decrease in the dephasing time for in- 
creasing Kq, a sharp dip around Kq = 2.4 is visible. In 
this region, J is suppressed by a factor of more than 20 
and hence the effective tunneling rate \Jce/h\ ^ 10 Hz, 
which for our experimental parameters is comparable to 
the on-site interaction energy U expressed in frequency 
units (we checked that the widths of the on-site wave- 
functions and hence U were independent of Kq by ana- 
lyzing the side-peaks in the interference pattern). This 
means that neighbouring lattice sites are effectively de- 
coupled and the local phases evolve independently due to 
interatomic collisions, leading to a dephasing of the ar- 
ray [ISi, ilTj . By increasing the dipole trap frequency (and 



hence U), we verified that the timescale for this dephas- 
ing decreases as expected. We also studied a re-phasing 
of the BEC when, after an initial dephasing at Kq = 2.4, 
the value of the shaking parameter was reduced below 
2.4. The time constant Treph of the subsequent rephasing 
of the condensate (mediated by inter-well tunneling and 
on-site collisions) increased with decreasing Joff (see the 
insert of Fig. 3 (a), where we compare Troph as a func- 
tion of Kq with the inverse of the generalized Josephson 

1 1/2 

frequency Wj^jsephson "^ff predicted by the two- well 
model EHoj)- 

Finally, we investigated the dependence of the dephas- 
ing time on the shaking frequency lu (see Fig. 3 (b)). In- 
terestingly, while the tunneling suppression as observed 
in-situ works even for two / J « 1, in order to maintain the 
phase coherence of the condensate, much larger shaking 
frequencies are needed. Indeed, for our system there ex- 
ists an optimum shaking frequency of fioj/ J « 30. 

In summary, we have measured the dynamical suppres- 
sion of tunneling of a BEC in strongly shaken optical 
lattices and found excellent agreement with theoretical 
predictions. Our results show that the tunneling sup- 
pression occurs in a phase-coherent way and can, there- 
fore, be used as a tool to control the tunnelling matrix 
element while leaving the on-site interaction energy un- 
changed (in contrast to the usual technique of increasing 
the lattice depth, which changes both) and without dis- 
turbing the condensate. This might ultimately lead to 
the possibility of controlling quantum phase transitions 
by strong driving of the lattice. In this context, it will 
be important to investigate the question of adiabatic- 
ity when dynamically changing the shaking parameter. 
Furthermore, our system also opens up other avenues of 
research such as the realization of exact dynamical local- 
ization using discontinuous shaking waveforms [1, [2l| or 
tunneling suppression in superlattices [22] . 
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